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For a quantum gas, being subject to continuous feedback of a macroscopic observable, the single-particle dy- 
namics is studied. Albeit feedback-induced particle correlations, it is shown that analytic solutions are obtained 
by formally extending the single-particle Hilbert space by an auxiliary degree of freedom. The particle's motion 
is then fed by colored noise, which effectively maps quantum-statistical correlations onto the single particle. 
Thus, the single particle in the continuously controlled gas follows a non-Markovian trajectory in phase-space. 



o 
o 

> 

o 
o 



> 
oo 
o 



o ■ 
^— > ■ 

G ■ 

3 : 



PACS numbers: 05.40.-a, 05.30.-d, 05.30.Jp 

As noted by Caves and Milburn 1 1], the continuous obser- 
vation of an object inevitably results in an increasing uncer- 
tainty of its momentum, which may cause it to eventually 
escape from the region of observation. An experiment thus 
necessarily requires a mechanism to keep the object at fixed 
location in the laboratory. Such constraining forces can be 
covered under the general notion of feedback control. 

Given limited observational capabilities, in most cases the 
macroscopic body is observed as a whole, its internal struc- 
ture remaining unresolved. Clearly the same holds also for a 
continuous feedback process, where measurements and con- 
ditioned unitary actions follow in continuous sequence. In 
consequence the issue may be raised, how the microscopic 
constituents of the object are affected by the feedback. 

It is well known that macroscopic systems with large num- 
bers of correlated constituents allow for exact solutions only 
for few model systems |2fl. The correlations, being respon- 
sible for the large degree of complexity, are usually gener- 
ated by interactions between the constituents. Since inter- 
nal and external degrees of freedom are usually considered as 
being decoupled, one might think that macroscopic observa- 
tions have no effect on the internal dynamics of a system 1 13]. 
However, strong and high-order correlations emerge also due 
to measurement or feedback control of a macroscopic observ- 
able n. 

In this Letter we show, that substantial effects on the 
particles of a continuously feedback-controlled quantum gas 
emerge. We prove, that, despite the above mentioned prob- 
lems of correlated systems, continuous feedback can be ex- 
actly solved for the single-particle dynamics. This dynamics 
contains effects due to particle correlations, that are formally 
introduced as colored noise feeding the single-particle motion. 

Continuous measurement 111 m and feedback have 
been studied and experimentally realized 1 7] in the past for 
various systems, such as single atoms or single harmonic sys- 
tems. Different from these treatments, we address here the 
case of a bosonic many-particle system and thus deal with the 
additional complexity due to particle correlations. 

Consider an ideal gas where the constituents of mass m are 
indistinguishable, bosonic particles. In order to keep at all 
times the center of the cloud of particles at a predefined tar- 
get position, a feedback loop is continuously applied to com- 
pensate for the motion of the center of mass of the gas. The 



feedback will act as a damping force and unavoidably also 
randomizes the motion of the center of mass. 

The quantum dynamics due to the continuous feedback pro- 
cess has been shown 1 1 ] to be governed by the master equation 
of quantum Brownian motion |8|,|9|], 
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Here H is the Hamiltonian of the free evolution of the many- 
particle system, which we specify as linear dynamics of an 
ideal gas, and the additional terms are due to the feedback. 

Continuous feedback is obtained as the limit of a discrete 
series of measurements and kicks of the center of mass X: 
In each discrete step X is measured with resolution ao and 
then a kick is applied leading to a shift — (oX, where X is 
the measured value. After that the free evolution takes place 
and the sequence is repeated. The continuous limit, where the 
average rate of feedbacks 7 — > 00, then requires that ao — > 00 
and ( — ► 0, such that a = ao/ yfy and C = Co7 remain constant. 
Thus the parameters a and £ in Eq. ([0 describe the strengths 
of measurement and subsequent kick. 

It should be emphasized that for the gas considered here, 
Eq. Q describes the dynamics of the many-particle density 
operator of the system g. Needless to say, that for the com- 
plete many-particle problem neither exact analytical nor nu- 
merical solutions exist. Nevertheless, it is straightforward to 
obtain solutions for collective variables, such as the center of 
mass X(t) or the total momentum P(t). In fact, these prop- 
erties will lack any features due to the many-particle aspect, 
but will be identical to that of a single-particle system [l|] of 
correspondingly larger mass. 

The center-of-mass rms deviation of a harmonically trapped 
gas with trap frequency ojq, for instance, will converge to the 
stationary value 
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Here AX = y/h/ (2McJo) is the ground-state width of the 
center of mass (total mass M) and 77 = (AXq) 2 / ((a 2 ) de- 
notes the ratio of spatial localization due to the potential over 
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that due to the feedback. The stationary value (0 is typically 
smaller than the initial value and thus indicates the gain in 
localization due to the feedback. 

Different from a consideration of the controlled variables 
themselves, in this Letter we address a more subtle question. 
It is the question, how the single particle in the gas is affected 
by the continuous feedback. Various important properties of 
the gas, that may be derived from the single-particle behavior, 
for example its density profile, justify this approach. 

The single-particles dynamics can in principle be deduced 
from Eq. Q, accompanied with the conceptual difficulties due 
to many-particle correlations. We proceed instead by mapping 
these correlations onto single-particle fluctuations, which will 
be shown to allow for analytical solutions. 

Let us first consider the feedback in more detail in a second- 
quantized picture, using the bosonic matter field <f>(x) with 
commutator relation [(j>(x), (jy{x')] =8(x—x'). The total mo- 
mentum of particles in the gas, as used in Eq. {0, then reads 



P = —ih f dxcbUx) 4>{x). 
J dx 
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It should be noted, that the operations of measurement and 
kick are typically implemented by use of external (e.g. op- 
tical) probe and control fields, that equally interact with all 
particles. Thus, the accessible properties of the gas are neces- 
sarily of extensive (i.e. additive) type. 

The center of mass of the gas, however, is an intensive 
quantity: It is the ratio of summed particle positions over 
the particle number. In conclusion this variable cannot be 
accessed in the above mentioned way. Instead, a priori in- 
formation on the particle number is required for rendering X 
quasi into an extensive quantity. Typically, the best informa- 
tion one can get on the particle number from some previous 
measurements may be its expectation value (N) . Using there- 
fore this average particle number, the truly accessible (quasi) 
center-of-mass reads 1 14] 
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Our goal is to gain information on the dynamics of a sin- 
gle particle in the gas, described by the time evolution of the 
reduced single-particle density matrix, 



p(x,x') = (ft{x')4>(x)). 



(5) 



Deriving the equation of motion of this density matrix from 
Eq. 0, an infinite hyrarchy of coupled equations for parti- 
cle correlations of increasing order would be obtained. This 
would not allow for a closed equation of motion for the single- 
particle density matrix (|5} alone. Of course, this is the central 
problem in many -particle physics, where only few model sys- 
tems provide analytical solutions 01 . In general one is forced 
to rely on approximations, such as the truncation of higher- 
order correlations. 

In our case, the many-particle correlations do not decay 
with increasing order, they represent the strong correlating ef- 
fects of the measurement of the collective variable X. Any 



outcome of the measurement can be realized by a vast num- 
ber of microstates compatible with the observed value X. Af- 
ter the measurement all these microstates become parts of the 
highly correlated many-particle superposition state. In addi- 
tion, the indistinguishability of bosonic particles leads to fur- 
ther quantum-statistical effects in this conditioned quantum 
state. 

There have been solutions to continuous-feedback prob- 
lems for many-particle systems, in effect however, consider- 
ing distinguishable particles II (ll . In this letter we deal with 
bosons, requiring thus a conceptually different approach. One 
should note, that due to the indistinguishability, Eq. (|5j is the 
expectable behavior of a representative particle in the gas and 
clearly not that of a specific (distinguished) one. 

Given the above mentioned problems, we proceed by ex- 
tending the Hilbert space beyond that of p(x, x'), such as to 
obtain closed equations of motion. A replacement for Eq. (0 
is searched for that, however, contains the complete informa- 
tion of p(x, x') and even more. A suitable candidate for such 
an approach is the following function 

W N (x,p;X,P) = (2ttH)- 4 J dx' J dX' J dP' J dS 
x exp{-i[x'p+PX'+XP'+(N-l)S] /h} 



It can be interpreted as the Wigner function of a joint two- 
body system, conditioned on N particles being in the gas. 
The first body of the system is given by a single particle with 
phase-space variables x and p. The second one is the center of 
mass of the "other" A— 1 particles with phase-space variables 
X and P. 

Tracing Eq. (jfjl over the macroscopic variables N, X, and 
P, the single-particle Wigner function, being equivalent to 
p(x, x'), is obtained. However, integrating over x, p one does 
not obtain the true center-of-mass Wigner function. For a 
large average particle number (N) 1, though, the annihi- 
lation of only a single particle by the field operators in Eq. (|6jl 
can be neglected. Then Eq. converges to the true two-body 
Wigner function of single particle and center-of-mass system. 

We can show that for Eq. ijfji a closed equation of mo- 
tion can be derived from Eq. Q. In consequence, also the 
single-particle state p(x,x') can be obtained by tracing over 
the macroscopic variables N, X, P. For particles of mass m, 
being bound in a harmonic potential of frequency ujo, for ex- 
ample, Wn(x,p; X, P) obeys the following Fokker-Planck 
equation 
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Here 6jv = (N-1)/(N), the total mass is M — m{N), and 
N plays the role of a parameter. This Fokker-Planck equation 
is of linear type with a positive semi-definite diffusion matrix. 
Thus a bound analytic solution for its Green function of Gaus- 
sian type can be found 1 11]. In consequence, the analytic so- 
lution for p(x, x') can be obtained in a rather straightforward 
way from the (analytic) solution of Eq. 0. 

For getting more insight into the dynamics described by 
Eq. 0, we now turn to the equivalent TV-parameterized 
stochastic differential equations, that read 
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dP N = 
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with £i, £ 2 being statistically independent Wiener processes. 

Besides the free evolution of the single particle and of 
the (quasi) center of mass, these equations also contain a 
feedback-induced coupling between the microscopic single- 
particle and the macroscopic center-of-mass degree of free- 
dom. Thus our method of solution, that incorporates all many- 
particle correlation effects in the single-particle dynamics, 
was obtained by allowing for a coupling with an auxiliary 
macroscopic degree of freedom. 

It is important to note, that since the drift matrix of Eqs 0- 
dl It is not of normal form, it is impossible to diagonalize this 
set of dynamical equations. Furthermore, both systems are 
fed by the same noise sources, which is apparent since both 
macroscopic and microscopic systems are subject to the same 
feedback process. 

When tracking the motion of both systems, i.e. microscopic 
and macroscopic, Eqs dSl- Jl 1> generate Markovian trajecto- 
ries in the four-dimensional phase space. The macroscopic 
system, however, being the center of mass of the N— 1 "other" 
particles, is not a physically accessible observable. Thus it 
appears natural to eliminate this auxiliary degree of freedom, 
to obtain stochastic differential equations for the microscopic 
trajectory alone. Let us perform this elimination for the case 
of large particle numbers: 

For (N) — > oo, due to its much larger inertia, the (quasi) 
center-of-mass system will no longer be affected by the mo- 
tion of the single particle. Thus in Eq. fl!0l > the term pro- 
portional to xn I (N) can be discarded, which decouples 
the macroscopic system. Consequently, the formal solution 
.Xjv(t) is then obtained as Ornstein-Uhlenbeck process. In- 
serting it into Eq. 0, one is finally left with equations of mo- 
tion for x 7v andpjv alone. 

The modifications of Eq. l|8} due to the elimination proce- 
dure are formally done by two replacements: Xn in Eq. 
becomes the deterministic part of the solution Xjy(t), which 
now acts as an external drive. And the Wiener increment 
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FIG. 1: Noise spectrum for qjv < l/v2 (ajv = 1/2) for the case 
satisfying Eq. J14I with r\ = V2 (solid), and for r\ = 1/ V2 (dashed). 



is replace by <i£jv(£), which is now a sum of a Wiener 
and an Ornstein-Uhlenbeck process: 

d£ N (t) = d£i(f) - 2T N dt ^dCi{t')coa[n N (t-t')} 



(12) 



Here Q 2 N =ujI-T 2 n and a at =T n /ujq, with T N = (Q N /2 be- 
ing the feedback-induced damping of the coherent oscillation 
in the trap potential. 

To characterize this effective noise source we consider the 
Fourier transform of the stationary correlation function 



S N (u)= lim dTe™ T £ N (t+T)( N (t), 



(13) 



which is a symmetric spectrum depending on N only via 
the parameter ajy. From Figs ^ and [2] it is observed that 
only in the high-frequency limit the white-noise background 
Sn(u>) = 1 is reached. At low frequencies several extrema 
appear depending on the values of and 77 (or instead Q. 

For weak feedback-induced damping, ctN < l/v2, the 
peaks at the damped center-of-mass oscillation frequencies 
±J7at can be resolved under the condition 



(14) 



see Fig.Q Such large values of r\ represent strong feedback- 
induced localization of the center of mass. For smaller values 
of 77 the two peaks overlap and only a single maximum at uj = 
remains (see dashed curve in Fig. 0. A remarkable effect, 
however, appears as minima are formed at frequencies larger 
than uiq. These minima represent noise reduction below the 
white-noise background and are generated by self-correlation 
of white-noise with its modulated version, as seen from the 
contributions to Eq. dl 2i . 
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FIG. 2: Noise spectrum for ajv > l/\/2 (?7 = 1/V2), for the oscil- 
latory case qjv = 0.8 (solid) and the over-damped case a at = 1.25 
(dashed). 



For a at > 1/V2, i-e. for large damping, the maxima at 
±Ojv can never be resolved, leaving the single maximum at 
u) = 0. Also in this case symmetric minima at noise levels 
below unity appear at frequencies slightly higher than the trap 
frequency ujq, cf. Fig. [2] 

This noise reduction may be understood as follows: Both 
the single particle and the (quasi) center of mass are subject 
to the same feedback-generated noise. With some time delay 
the center of mass transfers its noise to the single particle via 
the coupling, cf. Eqs (|8} and ( I10i . Thus the effective noise 
seen by the single particle at the frequencies of the minima 
is partially compensated due to destructive phase shifts of the 
two paths of the noise input. 

Since the noise feeding the single-particle coordinate is 
now colored, as seen from the previous discussion, the re- 
sulting stochastic trajectories in single-particle phase space 
will be non-Markovian. This feature is due to the feedback- 
generated strong many-particle correlations, that are now cast 
into an effective noise source for the single particle. Thus the 
coloredness of the noise Mil is a manifestation of generation 
of correlations between the particles. 

One may interpret this also as a feedback-mediated ef- 
fective interaction of the single particle with the cloud of 
surrounding particles, that bears a non-vanishing correlation 
time. The latter may be quantified by the inverse spectral 
width of Eq. Jl 3b . Thus the memory effects in the particle 
cloud can be made responsible for the non-Markovian trajec- 
tory of the single particle. 

Summarizing, we have shown that a specific many-particle 
problem, that is dominated by strong particle correlations can 
be analytically solved for the reduced single-particle dynam- 
ics. The solution was obtained by first extending the Hilbert 
space of a single particle by an auxiliary degree of freedom to 
obtain analytically solvable closed equations of motion. Elim- 
inating the auxiliary degree of freedom, it could be shown 
that the single particle follows a non-Markovian trajectory in 



phase space. Thus the many-particle correlations have been 
mapped onto a coloredness of the noise feeding the particle's 
motion. 

The specific problem, for which this solution was obtain is 
quite general, in that it represents typical control of a system 
on the macroscopic level. In conclusion even such a macro- 
scopic intervention creates correlations and a resulting non- 
Markovian behavior of single particles inside the controlled 
system. 

One may ask for a generalization of our approach: Is it gen- 
erally possible to map many -particle correlations onto specific 
features of noise sources feeding the motion of single parti- 
cles? Clearly in the presented case, the correlations were not 
generated by true particle interactions but by continuous feed- 
back. Thus an immediate application to problems of interact- 
ing particles seems not obvious. Considering however, that 
we can interpret the solved system equally well as quantum 
Brownian motion of scattering, and thus interacting, particles, 
a potential route to a more general methodology may come 
into sight. 
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